o 

(N 



The Integrated Relativistic Iron Line from Active Galactic Nuclei: 
Chasing the Spin Evolution of Supermassive Black Holes 

D. R. Ballantyne 

Center for Relativistic Astrophysics, School of Physics, Georgia Institute of Technology, 

Atlanta, GA 30332 



d ' david.ballantyne@physics.gatech.edu 

ABSTRACT 

m 

hM \ The spin of a supermassive black hole (SMBH) is directly related to the ra- 

r^ \ diative efficiency of accretion on to the hole, and therefore impacts the amount of 

5^1 fuel required for the black hole to reach a certain mass. Thus, a knowledge of the 

O I SMBH spin distribution and evolution is necessary to develop a comprehensive 

'^ \ theory of the growth of SMBHs and their impact on galaxy formation. Currently, 

' the only direct measurement of SMBH spin is through fitting the broad Fe Ka 

line in AGNs. The evolution of spins could be determined by fitting the broad 

K»- \ line in the integrated spectra of AGNs over different redshift intervals. The accu- 

^ I racy of these measurements will depend on the observed integrated line strength. 

^ ' Here, we present theoretical predictions of the integrated relativistic Fe Ka line 

strength as a function of redshift and AGN luminosity. The equivalent widths 

O ' of the integrated lines are much less than 300 eV. Searches for the integrated 

, ■ line will be easiest for unobscured AGNs with 2-10 keV luminosities between 

44 < \ogLx < 45. The total integrated line makes up less than 4% of the X-ray 

L^ background, but its shape is sensitive to the average SMBH spin. By following 

;_i ■ these recommendations, future International X-ray Observatory surveys of broad 

Fe Ka lines should be able to determine the spin evolution of SMBHs. 

Subject headings: galaxies: active — galaxies: nuclei — galaxies: Seyfert — 
surveys — X-rays: diffuse background — X-rays: galaxies 



1. Introduction 

In astrophysics, black holes can be described by only two parameters: their mass 
Mbh and their angular momentum J (parametrized by the dimensionless spin parameter 
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a = cJ/GM'^Yi where < a < 1). Over the last several decades, estimates of black 



hole masses have been obtained by measuring the gravitationa 
nearby objects such as binary stellar companions ( e.g., BoltonJ 



Ghez et al.l 



influence of the hole on 



1972 : ICowlev et al.l Il983 
20081 ). st ellar cusps (e.g. 



McClintock fc RemillardI |2006[ ). orbiting stars (e.g., 

Magorrian et al.lll998l ) or dense clouds of ionized gas (e.g., IPeterson et al.l 120041 ). However, 
measurements of black hole spin are much more difficult to obtain as the spin only sig- 
nificantly affects the shape of s pace-time out to ^ Q Vg {r-g = GMbh/c^ is a gravitational 
radius) from the event horizon (JBardeen et al.lll972l ). In particular, the value of the spin is 
directly related to th e radius of the innermost stable circular orbit (ISCO) around the hole 
(JBardeen et al.lll972l ). and therefore impacts the radiative efficiency of thin accretion disks. 
If the accretion flow is rotating in the same sense as the hole, then a larger a brings the 
ISCO closer to the event horizon, increasing the efficiency that an accretion flow can convert 
gravitational potential energy into radiation. A larger fraction of the rest-mass energy is 
radiated away in a highly efficient accretion disk, reducing the growth rate of the black hole 
mass. Th erefore, measuring th e distribution of spins w ill elucidate the role o f steady accre- 



mergers (e.g.. 



tion (e.g.. lVolonteri et al.ll2005l). chao tic accretion (e.g.. King fc Pringldl2006l ) and black hole 



Berti fc Volonteril 120081 ) in the build-up of supermassive black holes (SMBHs). 



X-ray variability and optical microlensing data show that the X-ray source in AGNs 



is within 10 Vg (e.g., Ilwasawa et al.l |2004J : iDai et al.l l2010l ). indicating that the accretion 



disk is likely illuminated by X-rays down to the ISCO. As spin clearly impacts the radius 
of the ISCO, the extent of the red wing o f the relativistically broa dened Fe Ka line is a 



direct measurement of black hole spin (e.g.. lReynolds &: Fabianll2008l ). In recent years, long 



observations by XMM-Newton and Suzaku have used this metho d to obtain SMBH spin 
estirnates for a small nurnber of bright, nearby Seyfert 1 galaxies (JBrenneman fc Reynolds 



20061 : iMiniutti et al.ll2009l : ISchmoU et al.ll2009l ). These type of spin measurements depend 



on precise spectral fitting to the broad Fe Ka line and continuum and are most easily 
performed on sources with strong Fe Ka lines. Thus, the larger collecting area provided 
by the future International X-ray Observatory {IXO) will be necessary in order to extend 
these types of spin measurements over a significant range of redshift. In the meantime, 
attempts have been made to measure the average Fe Ka line profile by stacking the X-ray 
data of the many AGNs detected in deep Chandra or XMM-Newton surveys over a wide 
range of redshift and AGN luminosity. If a significant Fe Ka signal can be detected in the 
integrated spectrum then, in principle, the average spin of the black holes contributing to 
the line could be determined. Again, the abili ty to make this measurern ent depends crucially 
on the strength of the integrated Fe Ka line. IStreblyanska et al.l (|2005[ ) averaged the XMM- 
Newton spectra of type 1 and type 2 AGNs detected in the Lockman Hole and found a strong 
(equivalent width [EW] ~ 500 eV) and broad Fe Ka line that indicated a non-zero spin for 
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th e average SMBHs in thei r sample. However, similar stacking analyses of different datasets 



bv iBrusa et al 



(l2005h a nd ICorral et al.l (120081 ) were not able to detect a broad component, 



with ICorral et al.l (|2008[ ) finding that the EW of the broad component must be < 400 eV. 
These attempts to measure the integrated broad Fe Ka line are technically challenging and 
are limited by the unknown intrinsic distribution of broad Fe Ka line strengths. 



Recently, iBallantynd (120101 ) described a method by which the average EW of the broad 
Fe Ka line could be calculated as a function of z and luminosity. That paper concen- 
trated on the distribution of likely EWs from individual sources, and did not consider the 
result of integrating lines over specific redshift int ervals, a step cruc ial to measuring the 
spin evolution of SMBHs. This letter builds on the IBallantynd ( 120101 ) method to calculate 
the integrated Fe Ka line over several different redshift and luminosity ranges. We also 
consider the dependence of AGN obscuration on the appearance of the integrated Fe Ka 
lines. These predictions show which regions in the redshift and AGN luminosity plane will 
have the largest average broad Fe Ka EW, and will be vital for guiding current and future 
attempts to map out the spin evolution of SMBHs by fitting the average Fe Ka line profile. 
The calculations are reviewed in the next section, and the results are presented in Sect. |31 
The following A-dominat ed cosmology is ass umed in this paper: H^ = 70 km s~^ Mpc~^, 
Qa = 0.7, and Qm = 0.3 JSpergel et al.ll2003h . 



2. Calculations 



The calculation of the mean Fe Ka EW as a function of z and A GN 2 -10 keV lum ir iosity , 
Lx, broadly follows the procedure described by IBallantynd ( l2010l ). The lUeda et al.l (120031 ) 
hard X-ray luminosity function ( HXLF) is com bined with the observed AGN black hole mass 
function (BHMF) at z = 0.15 (JNetzerl |2009| ) to calculate an Eddington ratio distribution 
at a given Lx- In order to predict the integrated Fe Ka line over several cosmo logically 
interest ing redshift ranges we assume that the BHMF evolves as (1 + 2;)^'^^ between z = 
and 5 (JLabita et al.ll2009l ). It is also assumed that the broad Fe Ka line is produced by 
a single strong reflection event from a dens e, geometrically thin accretion disk that ca n 
be approximated by a constant density slab ( jRoss fc Fabianlll993l : iBallantyne et al.ll200ll ). 
The EW of the iron line depends on the photon index of the illuminating power-law F, the 
relative abundance of iron Ap^, the ionization parameter of the disk ^, and the reflection 
fraction, R, which determines the relative strength of the reflected spectrum in the total 
observed spectrum ( Ballantyne et al.ll2002l ). The first three of these parameters are known 
to depend on the Eddin gton ratio of the AGN (JNetzer &: TrakhtenbrotI l2007l : llnoue et al. 
20071 : iRisaliti et al.l |2009|) . Refiection spectra are then computed for each Eddington ratio 
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contributing to an observed Lx- The Fe Ka EW is measured for each spectrum (assuming 
R = 0.5, 1, 2 or 4) and averaged over the Eddington ratio distribution to give a single value 
for that Lx- In order for the integrated line to reflect the changes in ^ and F, we keep track 
of the average energy at which the Fe Ka line reaches a maximum (-Emax) • These calculations 
are performed for 100 values of Lx between logLx = 41.5 and logLx = 48, and for 100 
values of z between z = and 5. 

These results can now be included in a X-ray background (XRB) synthesis calculation 
to predict the strength and profile of the integrated broad Fe Ka line over several differ- 
ent ranges_of_Juimnosity and r edshift. The synthesis model used here was last described 



by [Draper fc Ballantynd ( 120091 ). although we neglect the contribution of blazars for this 



application. The integrated rest-frame spectral intensity I{E) are computed by evaluating: 

r.r,. c r- f'^^^^T^ d^Lx.z) SEiLx,z)d^ 

^ ^"^1. Llt^ dlogLx (1 + znnUl + zf + fi^)V2^1°S^^c^^, (1) 



where d^{Lx,z)/d\ogLx is the lUeda et al.l (120031 ) HXLF, Se{Lx,z) is the absorbed rest 



frame spectrum of an AGN with intrinsic luminosity Lx at redshift z, and di is the luminosity 
distance to redshift z. At each {Lx, z) pair in the integral, a relativistically broadened Fe Ka 
line appropriate for a maximally-spinning (i.e., a ~ 1) SMBH is added to the spectrum Se 
such that it has the EW predicted by the reflection calculations described above. The 



relatiyistic profile for the line was calculated using the 'laor2' model in XSPEC12 (JLaor 



1991 



(cf. 



Arnaud 



Nandra et al 



1996) with a line emissivity of between 1.2 and 6 Vg and —3 from 6 to 400 Vg 



20071 ). The inclination angle of the line emitting material was taken to 
be 30 degrees, and the rest energy of the line is set to (-Emax) to account for changes in the 
average ionization state of the line as Lx and z are varied. Finally, to be consistent with 
the calculation of the Fe Ka lines, the photon index of Se{Lx, z) is derived from averaging 
F over the Eddington ratio distribution for a given {Lx,z) pair. This method translates 
to average photon indices varying from 1.4 to 2.5 depending on the value of z or Lx- In 
practice, this method of determining F for Se{Lx, z) has a negligible effect on the results. 

The integrated AGN spectra are then computed for five different redshift ranges: < 
-2<5, 0<z<l, 1<2;<2, 2<2<3 and 3 < 2; < 5. For each redshift range, five intervals 
in luminosity are considered: 41.5 < logLx < 43, 43 < logLx < 44, 44 < logLx < 45 and 
45 < logLx < 48. The EW of the Fe Ka line is measured by directly integrating the resulting 
spectra. The EW measurements are performed separately for type 1 AGNs (those with 
hydrogen column densities logA^n < 22) and type 2 AGNs (those with 22 < logA^n < 25), 
although the EW of the broad line in this model is independent of the obscuration toward 
the AGN. The impact of obscuration will be discussed at the end of Sect. 13. 1[ 
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Results 



3.1. Integrated EWs 



Figured] plots the measured EWs from the integrated spectra described in the previous 
section. The sohd hues show the results for R = 1 while the dashed and dotted lines assume 
R = 2 and 0.5, respectively. The colors denote the different redshift ranges over which the 
integration was performed. The measured EWs in all cases are significantly less than 400 eV, 
consistent with the upper-limit derived by ICorral et al.l (|2008[ ). Indeed, for R= 1, the EWs 
are between ~ 90 eV and ~ 220 eV. In fact, even if i? = 2 for every AGN, the EW of the 
integrated Fe Ka line does not break 320 eV. If the average R of AGNs is closer to 0.5 then 
the majority of the integrated broad lines have EWs< 120 eV. 



In the model of Fe Ka production described here and by iBallantynd (120101 ) the EW of 
the line is ultimately a function of the Eddington ratio of the A GN. If the SMBH is a ccreting 
too weakly, there is no optically thick disk for reflection to occur (JNarayan fc Yilll995[ ). In the 
opposite extreme, if the accretion ratio is very close to the Eddingtoi i limit there is also littl e 
observed reflection signature because the disk is very highly ionized (JBallantyne et al.ll2002l ). 
The strongest Fe Ka line arises from recombinati on onto He-like iron, bu t this occurs over 
a relatively narrow range of ionization parameter (iBallantvne et al.l 120021 ). and, therefore, a 
correspondingly small range of Eddington ratios ( lB allantynd |2 1 Ol ) . This interplay between 
Eddington ratios and Fe Ka EWs is reflected in the redshift distributions of the integrated 
EW seen in Figure [U For all AGNs with log Lx < 45 the largest integrated EW is found 
from those with z < 1. As z increases, the average mass of active SMBHs increases, which, 
at a constant Lx, corresponding to a decreasing Eddington ratio. Therefore, the integrated 
Fe Ka EW drops as z increases because a larger fraction of AGNs at these luminosities 
are accreting too weakly to produce a strong disk line. Integrating over the entire redshift 
range (black lines in Fig. [T]) results in a weighted average of this behavior. In contrast, the 
Fe Ka EWs of AGNs with logLx > 45 are largest for 1 < z < 2 and smallest for z < 1. The 
luminosities of these AGNs are large enough that ionized Fe Ka lines contribute significantly 
at z > 1 and boost the integrated EW. In this case, the Eddington ratios at lower z are, on 
average, too large and the resulting disk lines are weaker. At larger redshifts, the increasing 
SMBH mass reduces the average Eddington ratio which results in strong neutral Fe Ka 
emission. This discussion illustrates that, for a fixed luminosity range, the integrated line 
probes smaller Eddington ratios as the redshift increases. Elucidating the spin evolution of 
the SMBH population over a fixed interval of mass will require measuring the integrated 
Fe Ka line over ranges of progressively higher luminosity as the redshift increases. 



At face value, the results of Fig. [T] indicate that experiments using the broad Fe Ka 
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line to measure the spin evolution of SMBHs should be targeted toward integrating the 
spectra of AGNs with logLx > 45, as these sources are the ones predicted to have the 
strongest integrat ed EWs. The diff iculty with this strategy is that such AGNs are rare at 



all redshifts (e.g., lUeda et al.ll2003[ ). and thus the compilation of suitable samples will be 
extremely problematic even with the enhanced sensitivity of IXO. Therefore, we conclude 
that AGNs with 44 < log Lx < 45 will provide the best sample to search for the integrated 
broad Fe Ka lin e. These AGNs are close to the knee of the HXLF at all redshifts (e.g.. 



Ueda et al.ll2003l ) and will therefore be common enough to provide a useful sample at several 
different redshift bins. The broad lines from these AGNs will show a mixture of of neutral 
and ionized reflection with (-Emax) ~ 6.5-6.6 keV. 

Examples of the integrated R = 1 spectra in the 44 < log Lx < 45 range are shown as 
solid lines in Figure [2] with the colors denoting the redshift ranges as in fig. [T] These spectra 
were calculated assuming that the Fe Ka line was being emitted all the way down to the 
ISCO of a maximally spinning black hole. The dashed lines show the integrated spectra if 
the ISCO is at 6 Vg, as would be the case for a non-spinning black hole. As expected, this 
scenario causes the integrated line to be slightly less broad and more peaked. This small 
difference in the line profile should be measurable if a large enough sample of AGNs could 
be included in the integral. This result nicely illustrates the importance of targeting AGN 
in the luminosity range that will produce the most intense integrated broad Fe Ka line. 

The integrated spectra shown in figure [2] and the associated EWs plotted in fig. [T] are 
compiled from only type 1 AGNs. The dotted blue line in Fig. [2] plots the integrated type 
2 spectrum with 44 < logLx < 45 and z < 1. Our calculations assume that there is no 
difference in the broad Fe Ka EW between the type 1 and 2 AGNs. However, a visual 
inspection of the type 2 spectrum in fig. [2] seems to indicate a much more intense broad 
Fe Ka line. Indeed, performing the same integration of the type 2 line gives an EW of 
460 eV as compared to 180 eV for the type 1 AGN. The strong absorption in the average 
type 2 spectrum causes the continuum to curve downwards at energies < 5 eV and enhances 
the Fe K edge at 7.1 keV. These effects conspire to artificially enhance the EW of the broad 
Fe Ka line. Therefore, unless the underlying spectrum and absorption distribution is well 
known, the search for broad Fe Ka lines in integrated spectra should only be performed on 
unobscured, type 1 AGNs. 



3.2. The Broad Line in the X-ray Background 

Performing the integral of eqn. [T]over all z and logLx, but now redshifting Se into the 
observed frame results in a prediction of the contribution of the Fe Ka line to the entire XRB. 
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Other authors have made similar predictions ( IGiUi et aLlll999l : iGandhi fc Fabian! |2003| ). but 
have had to make very simple assumptions on the Fe Ka strength and shape. Figure [3] shows 
the first prediction of the Fe Ka contribution to the XRB that arises from a physical model 
of the Fe Ka line distribution. This figure plots the ratio of the XRB spectrum including 
the broad line to one which does not include the broad line. Predictions are shown for R = 1 
(solid line), 2 (short-dashed line) and 4 (dot-dashed line). The long-dashed line plots the 
ratio for a model where the broad line has a constant EW of 100 eV and results from neutral 
reflection. In all of these cases the line is broadening using the 'laor2' model as described in 
Sect. [21 except for the dotted line which plots the R = 1 model with all lines arising from 
6 Vg (also using the 'laor' model). 

Figure [3] shows that the maximum contribution of the broad Fe Ka lin e to the XRB 



is 6%, but, since inten se reflection cannot be a common occurrence in AGNs (JGandhi et al. 
20071 : lBallantyndl2010l ) , the contribution is more rea listically rnuch les s tha n 4%. These value s 



are in agreement with the previous predictions of iGilli et al.l (119991 ) and iGandhi fc Fabian 
(120031 ). The contribution of integrating multiple observed- frame spectra and relativistic 
blurring smears the Fe Ka contribution over 6 keV in energy, but there is a slight difference 
between the profiles with different spin distributions. Therefore, given a very high qual- 
ity spectrum of the X-ray background, the total integrated Fe Ka profile will be able to 
determine the average spin of SMBHs over cosmic time. 



4. Conclusions 

Measurements of the spin distribution of SMBHs will have profound consequences on 
our understanding of galaxy formation and evolution. The only direct way to measure SMBH 
spins is via the broad Fe Ka line in AGNs, and attempts have been made to measure the 
Fe Ka profile in the spectrum of AGNs that have been integrated over a range of z and 
logL^. The accuracy and precision of these measurements will be maximized if the spectral 
fitting can be performed on Fe Ka lines with the largest possible EW. Therefore, this letter 
presented predictions of the integrated relativistic Fe Ka EW as functions of z and logL^. 
Fig. [1] showed that the integrated broad Fe Ka EW will be significantly less than 300 eV for 
all ranges of z and logLx. In fact, if i? = 1 for most AGNs, the integrated broad component 
will have an EW between 100 and 200 eV. Taking into account the relative space density 
of AGNs, it is recommended that the search for integrated broad Fe Ka lines be confined 
to AGNs with luminosities in the range 44 < logLx < 45. It is further recommended that 
these experiments be limited to type 1 AGNs. 



Stepping back from the quantitative results, this letter has shown that measuring the 



spin evolution of SMBHs can be done with the broad Fe Ka hne. In fact, IXO will be 
perfectly placed to obtain this goal, as its sensitivity will allow it to obtain the high-quality 
spectra necessary to measure the broad Fe Ka lines. IXO spectral surveys of target lists 
provided by the Chandra deep fields will be a very powerful tool to trace the spin of SMBHs 
up to and beyond z = 1. 

The referee is acknowledged for a detailed and useful report that greatly improved the 
paper. 
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Fig. 1. — Fe Ka EWs measured from AGN spectra integrated over different ranges of z and 
Lx (Eqn. [1]). The colors denote the different redshift ranges the integrals were performed 
over with black indicating the result from the entire 2; = to 5 interval. The solid lines show 
the results for a reflection fraction of unity (i.e., R = 1), while the dashed and dotted lines 
plot the EWs if all AGNs have R = 2 01 0.5, respectively. 
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Fig. 2. — Examples of the integrated R = \ spectra from which the Fe Ka EWs were 
measured. The colors correspond to the same redshift ranges as in Fig. [1] The solid lines 
show the spectra for the 44 < logLx < 45 bin assuming the emission line extends down 
to the ISCO of a rapidly spinning black hole. The dashed lines plot the same spectra, but 
now assuming that the ISCO is 6 r^ for all AGNs. Both the solid and dashed lines are the 
integrated spectra including only type 1 AGNs (i.e., those with logA^n < 22). The dotted 
blue line plots the integrated type 2 spectrum for AGNs with 44 < logLx < 45 and z < 1. 
The lines in this spectrum also assume emission down the ISCO of a rapidly spinning black 
hole. In our model the EW of the broad Fe Ka line is independent of nuclear obscuration. 
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Fig. 3. — The ratio of a model XRB spectrum that includes the broad Fe Ka line to one 
that does not include the broad line. The ratios shown include a model with R = \ (solid 
line), R = 1 (short-dashed line) and R = 4 (dot-dashed line). Also shown is the ratio from 
a model where each AGN has a neutral broad Fe Ka with a constant EW of 100 eV (long- 
dashed line). Finally, the dotted line corresponds to a i? = 1 model where all the SMBHs 
are non-spinning and the lines extend to only 6 Vg. 



